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From the observation of a millisecond pulsar in SN 1987A, the following implications 
are obtained. 1) The pulsar spindown in SN 1987A is caused by radiating gravitational 
waves rather than by magnetic dipole radiation and/or relativistic pulsar winds. 2) A mildly 
deformed shock wave would be formed at the core-collapse and explosion in SN 1987A, which 
is consistent with the conclusion given in Nagataki (2000). 3) The gravitational waves from 
the pulsar should be detected in several years using a Fabry-Perot-Michelson interferometer 
as the gravitational detector, such as LIGO and TAMA. 4) The neutrino oscillation model 
is not a promising model for the explanation of the kick velocity of the pulsar in SN 1987A. 
The hydrodynamical instability model is more favored. 

§1. Introduction 

Since there is little informations concerning the pulsar in the remnant of SN 
1987A, its properties, such as angular velocity of rotation, strength of magnetic fields, 
and total baryon mass, have been treated as free parameters or output parameters. 
EP □) However, Middleditch et al. (2000) reported the discovery of an optical pulsar 
whose frequency is 467.5 Hz and spindown rate is (2-3) xlO -10 Hz s _1 .EP Since some 
free parameters appearing in previous papers are constrained by this discovery, we 
consider its implications in this paper. In section we show that the spindown 
is caused by radiating gravitational waves rather than by magnetic dipole radiation 
and/or relativistic pulsar winds. We also determine constraints on the strength of the 
magnetic field of the pulsar. In section ||, we discuss the effects of the proto-neutron 
star's angular momentum on the dynamics of the core collapse of the progenitor of SN 
1987A. The amplitude of the gravitational waves from the pulsar and its detectability 
are discussed in section H. Implications of the kick velocity of the newly-born pulsar 
are presented in section || 

§2. Origin of the pulsar spindown 



Middleditch et al. (2000) reported that the spindowns (2-3xl0~ 10 Hz s _1 ) of 
the 2.14 ms pulsations should be caused by radiating gravitational waves. This 
is because the relation between the spindown rate and its modulation period can 
be explained at the same time by adding the non-axisymmetric component of the 
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moment of inertia (81) to a spherical neutron star whose moment of inertia is /. 
It is true that this conclusion is curious, because the spindown of a normal pulsar 
is believed to be caused by magnetic dipole radiation^ and/or relativistic pulsar 
winds. LP However, we want to emphasize that their conclusion is supported by the 
recent UVOIR bolometric light curve. We can easily calculate the decreasing rate of 
the rotational kinetic energy of the pulsar as 

alt di ' y ' 

where Q is the angular velocity of the pulsar. Assuming that the pulsar is spherical 
and has a constant density, the moment of inertia of the pulsar can be expressed as 

[gem 2 ]. (2-2) 



1.4 M Q ) UOkrn,/ 
Thus, using the observation of Q and f2, we can estimate E as 

E = -(4-6)xW™(-^)(-J-) 2 [ergs- 1 ], (2-3) 



.1.4 M & J \ W km, 

which is much larger than the UVOIR bolometric luminosisy, (1-2) xlO 36 erg s 1 . © 
This discussion strongly supports that the pulsar spindown is caused by radiat- 
ing gravitational waves. Otherwise, the supernova remnant would become much 
brighter. 

We must also check whether the remnant is bright in other wavelengths, such as 
radio, X-ray, and gamma-ray. If the brightness of the remnant in these frequencies is 
not too large, we can confirm more strongly our hypothesis that the pulsar spindown 
is caused by radiating gravitational waves. 

It is reported that the radio emission spectrum is well fitted as& 

S~w- l5 ( [ergs" 1 cm" 2 GHz- 1 ]. (2-4) 

V 1 GHz / 

In obtaining this expression, Gaesnler et al.H* used data at frequencies of 1.4, 2.4, 
4.8, and 8.6 GHz. Therefore, when we assume that this power-law fitting holds at 
all radio frequencies, we can estimate its luminosity as 

W ~ 3 X 10 32 ( ^-flGge ( , (2-5) 



50 kpc 

where D and v are the distance from the Earth to the remnant and the radio fre- 
quency, respectively. Assuming that the distance is 50 kpc, IHP we find that the 
luminosity in the radio band is much smaller than the rate at which the rotational 
energy decreases. In fact, unless f m in is as small as 10 -10 Hz, the luminosity in the 
radio band is not comparable to the rate at which the rotational energy decreases. 
Moreover, it is generally believed that the radio emission does not come from the 
pulsar but from the synchrotron emission of electrons that is generated when the 
shock encounters circumstellar matter. |) 
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An upper limit for X-rays of 2.3xl0 34 ergs -1 (0.5 - 2 keV) has been placed 
by the Chandra observations of the remnant. c3 They also discussed that this low 
upper limit is not surprising in view of calculations showing that debris should still 
be opaque to soft X-rays. Hence we can conclude that the rapid rate at which the 
rotational energy decreases cannot be explained by the emission of the soft X-rays. 

The remnant is thought to be transparent to hard X-rays and gamma-rays. In 
fact, Fransson and Chevalier t3 reported that the energy corresponding to unity of 
the absorption optical depth of the ejecta can be well represented by the formula 

/ M \ °- 36 / V \ -°- 72 ( t \ ~ - 72 

E(t = \)=%\[^-\ ( £ =-) ( — ) [keV],(2-6) 

1 ' \IM Q J v 2500 km s-V \lyr) 

where M c , V c , and t are the mass inside the O/He interface, the expansion velocity 
of the core, and the time from the explosion, respectively. When we adopt M c = 
3.7M ,§) V c = 2500 km s _1 ,B and t = 5 yr, we obtain E(t = 1) = 40 keV. Thus the 
situation here is different from that for soft X-rays; that is, the remnant is thought to 
be transparent to hard X-rays and gamma-rays. As for the data at these frequencies, 
the upper limit of the spectrum is rather rough, and published data are not new. 
The gamma-ray continuum from April 4, 1989 can be fit ast 



. Q3D 



= 1.6 x 10~ 5 [photons cm" 2 s" 1 keV" 1 ] (2-7) 

dE VIOOkeV/ 1 

for the energy range 50-800 keV. The total energy flux can be obtained as 



where E max is the maximum energy of the gamma-ray photons. It is generally 
believed that these gamma-rays come from radioactive nuclei such as 56 Co and Co. 
t3 Moreover, the Crab nebula, whose energy source is the central pulsar, is brightest 
in the X-ray band. Thus it is difficult to think that the remnant in SN 1987A could 
be brightest in the gamma-ray band and its luminosity could be comparable to the 
rate at which the rotational energy decreases. In the case that the total gamma-ray 
luminosity is found by future observations to be as large as the rate at which the 
rotational energy decreases, we will have to consider the serious problem of producing 
gamma-rays only. That is, we have to face the difficult problem of determining a 
mechanism that produces only gamma-rays and no photons in other energy bands 
from the pulsar. 

For the reasons mentioned above, we think that it is difficult to explain the 
observed spindown with the magnetic dipole model and/or relativistic pulsar wind 
model. However, it is necessary to fix a strict upper limit of the present gamma-ray 
flux from SN 1987A in order to conclude with certainty that the pulsar spindown is 
not caused by radiating photons and/or ejecting relativistic particles but by gravi- 
tational waves. 

We can also give un upper limit for the strength of the magnetic field of the 
pulsar. From the magnetic dipole model and/or the relativistic pulsar wind model, 
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the rate at which the energy decreases can be written 



I) 



E= — , (2-9) 



c 3 

where B p is the strength of the magnetic field at the magnetic pole of the star. From 
Eqs. (2-3) and (2-9), we can derive the upper limit for B p as 

Since the strength of the magnetic field is so weak, we believe that a neutron star 
without hot spots will be found in X-ray bands in the near future. The lack of 
hot spots implies that the surface temperature of such a neutron star will be ap- 
proximately uniform. Such a neutron star should be found when the optical depth 
becomes sufficiently low. We also add a comment on the weakness of the magnetic 
field of the pulsar in SN 1987A. Since it is apparently weaker than the typical one, it 
seems to be suggested that the strength of the magnetic field of a newly born pulsar 
evolves as a function of time. Thus we may be able to observe in the near future the 
magnetic field of a pulsar growing stronger when we can observe the pulsar activity 
directly. 

Finally, we consider the possibility of radiating gravitational waves to a suffi- 
cient extent to explain the observed pulsar spindown. Middleditch et al. (2000) con- 
cluded that the required non-axisymmetric oblateness (e) is ~ 10~ 6 , since a slightly 
deformed, homogeneous ellipsoidal pulsar with moment of inertia I and ellipticity e 
radiates energy in the form of gravitational 

£cw = -f ^/W. (2-11) 

o c 

Here e is defined as 

<=<£i^ < 2 ' 12 > 

where a and b are the equatorial semiaxes. 

We should discuss where and why a non-axisymmetric component of the moment 
of inertia is realized in a neutron star. Here we have to note that the average 
density of a pulsar is about 5xl0 14 g cm -3 . Thus it is meaningless to consider a 
'mountain' on the surface of a neutron star, where the density is about 10 9 g cm -3 
and the density scale height is only ~ 1cm.© This is because the contribution of 
the mountain on the surface of the neutron star is too little to explain the non- 
axisymmetric component of the moment of inertia. Rather, we should consider 
density fluctuation in the inner crust, where the typical density is sufficiently high 
and the contribution to the moment of inertia is not negligible. In particular, Lorenz 
et al.l!3 reported that there may be a nuclear 'pasta' at the innermost region of the 
inner crust. HIP In this nuclear pasta region, we can easily guess that non-uniform 
crystallization due to the rapid cooling of the newly born neutron star will result in 
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such a non-axisymmetric component of the moment of inertia. It is a very important 
task to estimate the nucleation rate and the growth rate of such a crystal in the 
pasta region. This should give information on the non-axisymmetric component of 
the moment of inertia in neutron stars. 

§3. Implications for the jet-like explosion in SN 1987A 



Effects of rotation on the dynamics of collapse-driven supernovae have been 
investigated in many works. 0) ' l§ ' li^P ' EE* However, the initial rotational energy has 
been given parametrically in these works, because we have little information on 
it. Because we now have information on the rotational energy of the newly born 
pulsar in SN 1987A, we can now carry out further analysis. First, we estimate the 
initial period of the pulsation. From Eqs. (2T) and (2-11), the initial period can be 
estimated as 

T i = T (l-4^), (3-1) 

where T , t, and Tqw are the present period (~ 2.14 ms), present time ~ 5 yr, 
and -J? /J? , respectively. From Eq. (3-1), the initial period of pulsation can be 
estimated as (1.9-2.0)ms. Even if we assume that the pulsar spindown is caused 
by magnetic dipole radiation and/or relativistic pulsar winds, the estimated initial 
period changes little. Now we can estimate the ratio of the rotational energy relative 
to the gravitational binding energy at the moment of the core-collapse (T/|W^|i n it). 
It is estimated as 

T/\WU= ^ (-) (3-2) 
, / M \ /1000 km\ 9 , . 

where q = Jc/GM 2 = IQc/GM 2 is the dimensionless angular momentum. Since 
the value for q can be estimated as 

T/\W\ init for the progenitor of SN 1987A can be estimated to be 1.7 xlO" 4 . 

We wish to stress the fact that this estimated value is smaller than the values 
assumed in the study of Yamada and Sato (1994), in which an extremely deformed 
shock wave is formed. Thus, it can be easily guessed that a mildly deformed shock 
wave was formed in the core of SN 1987A. This is consistent with the conclusion 
reached in Nagataki (2000). 

We can also estimate the ellipticity (e) of the proto-neutron star in SN 1987A 
from the rotational energy. Since the relation between e and T/|W| can be written 
as§> 



\W\ 




e(l-e 2 ) 1 / 2 ' 
sin -1 e 



(3-5) 
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e can be estimated as ~ 0.25. Here we have assumed that the mass and radius of the 
proto-neutron star are 1.4M© and 20 km, respectively. This means that the ratio 
of the semimajor axis relative to the semiminor axis of the proto-neutron star is ~ 
1.03. It should be noted that this value is smaller than that assumed in the work of 
Shimizu et al. (1994), in which an extremely deformed shock wave is formed due to 
the effects of asymmetric neutrino heating from the deformed neutrinosphere. This 
discussion given here also supports the conclusion reached in Nagataki (2000) that 
a mildly deformed shock wave is required in order to realize the appropriate matter 
mixing and explosive nucleosynthesis in SN 1987A. 

It is a very important task to perform numerical simulations in which the effects 
of rotation and neutrino heating are included in order to make an appropriate model 
for SN 1987A in which a mildly deformed shock wave and a rotating neutron star 
with a period of 2 ms are formed. Such a model will help us to understand more 
clearly the system SN 1987A and the roles of rotation and asymmetric neutrino 
heating on the dynamics of collapse-driven supernovae. 

§4. Gravitational waves from the pulsar 



We can estimate the amplitude of the gravitational waves from the pulsar in SN 
1987A. The energy release rate due to the gravitational waves can be written 

E= T ^0 2 (h) 2 x4ttD 2 , (4-1) 

where {h) is the average dimensionless amplitude of the gravitational waves at the 
distance D from the pulsar. From Eqs. (2-11) and (4T), (h) can be estimated as 

W-S.l^Ie^ 2 (4-2) 



~ 4.7 x 10" 26 



1.1 x 10 45 gcmV VIO- 6 / \2936 rads-V V D 



I \ ( e \ / Q V /50kpc 



-^4,3) 



Thus, the required time to detect the gravitational wave from the pulsar using the 
Fabry-Perot-Michelson interferometer as the gravitational detector is 

AT = 4.2 ( =] fyrl. (4-4) 

^3x io- 22 [i/VHz]y u J v ' 

where h! is the sensitivity of the detector at 2x467.5 = 935 Hz, which has units 
of l/\/Hz. We can find that detection of the gravitational wave from the pulsar is 
possible_within a reasonable time when gravitational detectors such as LIGOB1P and 
TAMAE3), whose sensitivities are of order /i'~10 -22 Hz -1 / 2 , are running. 



§5. Implications on the kick velocity of a newly-born pulsar 
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It is a well known fact that pulsars in our galaxy have velocities much greater 
than those of ordinary stars. c3> It is reported that their transverse speeds range 
from to ~ 1500 km s _1 and their mean three-dimensional speeds are ~450 ± 90 
km s 1 . IH^P There are many theoretical models to explain the pulsar kick. According 
to one, a neutron star in a binary system can escanairom the system with rapid speed 
due to a supernova explosion of the nascent star.E^ There are also many models in 
which effects of asymmetric supernova explosions are taken into consideration. HEZP 
For example, it is reported that neutrino oscillations, biased by the magnetic field, 
alter the shape of the neutrino sphere in a cooling proto-neutron star and are the 
origin of the kick velocity of the pulsar. 0) In another work, Burrows and Hayes 
(1996) pointed out the possibility that hydrodynamical instabilities may be the origin 
of the pulsar kick. However, there are too few observations to determine which model 
is the most promising one. 

It is suggested that the pulsar in SN 1S87A alsohas a kick velocity and is moving 



toward the south region of the remnant. E§) • r3 > [§>B As discussed in section ^, it is 
suggested that the strength of magnetic fields on the surface of the pulsar in SN 
1987A is very weak. Therefore it is concluded that a neutrino oscillation model like 
that of Kusenko and Segre (1996), which requires magnetic fields of order 10 14 G, is 
not promising. This is the first investigation with the purpose of selecting the best 
model of the kick velocity using observational data. We will be able to give further 
discussion when we obtain more precise data on the pulsar in SN 1987A. We hope 
there will be further observations of this pulsar at many frequencies of photon and 
gravitational waves in the near future so that we can continue our investigation and 
confirm the report of the discovery presented by Middleditch et al. (2000). 
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